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The Landau- Ginzburg- Wilson hamiltonian is studied for different values of the pa- 
rameter A which multiplies the quartic term (it turns out that this is equivalent to consider 
different values of the coherence length £ in units of the lattice spacing a). It is observed 
that amplitude fluctuations can change dramatically the nature of the phase transition: 
for small values of A (£/a > 0.7) , instead of the smooth Kosterlitz-Thouless transition 
there is a first order transition with a discontinuous jump in the vortex density v and a 
larger non-universal drop in the helicity modulus. In particular, for A sufficiently small 
(f/a = 1) , the density of bound pairs of vortex-antivortex below T c is so low that, v 
drops to zero almost for all temperature T < Tc. 

Vortices play a central role in explaining the phase diagram and properties of superfluid systems 
both neutral and charged. The discovery of high-temperature superconductors has boosted the interest 
in the dynamics of the vortex lines in the mixed state and opened a new area which regards the physics 
of vortices as a new state of matter |5J]. 

Superfluid films and Josephson junction arrays in two dimensions are often described by XY-type 
models in which the unique microscopic degrees of freedom are phases. Thanks to the work of Berezinskii 
Q and Kosterlitz and Thouless J3[ in the early 70's we have a fair understanding of the standard XY 
model in two dimensions. In their theory a thermodynamic phase transition, the Kosterlitz-Thouless (KT) 
transition is driven by the unbinding vortices ( singular phase configurations ) at a temperature Tkt- In 
refs. |Q- Q it was pointed out that the KT phase transition might also apply to thin-film superconductors 
and some experimental evidence was discussed. In particular, the analysis of the 2 dimensional flux line 
lattice (FLL) melting transition of ref. || is based on a KT-type theory. More recently, it was suggested 
that some features of the KT theory might be present in under-doped superconducting cuprates tt. 

However, the nature of the phase transition driven by vortices in 2D still remains under discussion. 
By means of Monte Carlo simulations of the XY model with a modified nearest neighbor interaction 
it was shown that, depending on the value of an additional parameter, continuous as well first-order 
transitions take place M]. The existence of both kinds of phase transitions is in accordance with 
the richer structure of the 2D Coulomb gas found by Minnhagen and Wallin (ll) using self-consistent 
renormalization group equations and with the tendency towards first-order transition which develops in 
the case of a strong disorder coupling constant fll2f . 

Our approach to the nature of the phase transition driven by vortices in 2D superflu- 
ids/superconductors is a different one based on the effect of amplitude fluctuations. Let us start by 
noting that right at the ( expected vortex-pair unbinding ) transition the amplitude fluctuations cannot 
be considered as weak anymore and thus may affect the critical behavior |l3f| . Therefore, instead of the 
XY model -with fixed amplitude fields- it is worth investigating the effects of vortices when amplitude fluc- 
tuations are not neglected. Hence, in this letter we analyze the more general Landau- Ginzburg- Wilson 
(LGW) lattice Hamiltonian, in terms of a complex field ip = |^>|expi# (IV'I ^ constant) and two 
parameters K and A p5| : 

0H = \^\\^+^os{e x+ail - e x ) + ]T[A(|^| 2 - 1) 2 + m 2 ], (1) 
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where f3 — 1/T, a is the lattice spacing, x denotes the lattice sites and the index fi — 1,2. We will 



show that the nature of the phase transition of this model diverges dramatically from the KT when the 
parameter A is chosen sufficiently small -which in fact is equivalent to take the coherence length £ ~ a- 
and that this is connected with the appearance of a sharp jump in the number of vortices. 

A straightforward discretization of the continuum Landau-Ginzburg Hamiltonian produces the ex- 
pression Q: 

[3H = (3a 2 Y,d2 t^Ua» ~ r x f/a 2 + r|^| 2 + u|^| 4 ], (2) 
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where the superscript c in ip denotes the ordinary parameterization in the continuum theory m is the 

effective mass of the carriers and the coefficients r and u are analytic functions of the temperature, with 

_ i 

u > for stability. Introducing a dimensionless order parameter: tp x — (^r) 2 ipx and writing 

0H = ^a 2 £5>x + a„ " ^)V« 2 + 2^(1 - I4I 2 ) 2 ], (3) 
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where = ^muP and £ is the coherence length given by £ 2 = 2 m\ r \ ■ Parameterizations (|l]) and (||) are 
connected by the relations: 
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In the limit of A = oo the radial degree of freedom is frozen and this model -sometimes said to 
describe soft spins with non fixed amplitude- becomes the XY model which is said to describe hard spins 
with fixed amplitude. A more interesting and less well studied limit is just the opposite i.e. small values 
of the A parameter. By (Q) small values of A correspond to a large £ in units of a i.e. A and l/£ 2 are the 
self- interaction coefficients respectively in (|l|) and (||) regulating amplitude fluctuations. 

We have simulated the Hamiltonian (pi) using a Monte Carlo algorithm. The calculations were 
performed on square L x L lattices with periodic boundary conditions (PBC). In order to increase the 
speed of the simulation we have discretized the 0(2) global symmetry group to a Z(N) and compared the 
results with previous runs carried out with the full 0(2) group in relatively small lattices. For the case of 
Z(60) we found no appreciable differences. Lattices with L = 10, 20, 24, 32, 40 and (in some cases) 64 were 
used. For L = 10, 20, 24 we thermalized with, usually, 20.000-40.000 sweeps and averaged over another 
60.000-100.000 sweeps. For L = 32,40 and 64 larger runs were performed, typically 50.000 sweeps were 
discarded for equilibration and averaged over 200.000 sweeps. We also performed some more extensive 
runs near K c for small values of A. The following quantities were measured: i) The vortex density v. The 
standard procedure to calculate the vorticity on each plaquette is by considering the quantity 

m = - e 2 ] 2v + \e 2 - 9 3 } 2 „ + [e 3 - e 4 ] 2v + \e A - 6i] 2n ), (5) 

where [a]2ir stands for a modulo 2tt: [a^K = a + 2nn, with n an integer such that a + 2un 6 (— tt, n], 
hence m = ni 2 + n 2 $ + 7134 + 7141. If m / 0, there exists a vortex which is assigned to the object dual 
to the given plaquette. Hence in the case d = 2, *m, the dual of to, is assigned to the center of the 
original plaquette p. The vortex "charge" *m can take three values: 0, ±1 (the value ±2 has a negligible 
probability), v defined as: 
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serves as a measure of the vortex density, ii) The energy density e =< H > / L 2 and the specific heat 
C v , which were computed to measure the order of the phase transition, iii) The helicity modulus T p6| 
which measures the phase-stiffness. For a spin system with PBC the helicity modulus measures the cost 
in free energy of imposing a "twist" equal to LS in the phase between two opposite boundaries of the 
system. T is obtained in general as a second order derivative of the free energy with respect to 5 -which 
can be regarded as a uniform statistical vector potential- evaluated for S — ► 0. In such a way one gets the 
following expression ]l7j , which generalizes the one introduced in ref. to an order parameter with 
amplitude as well as phase variations: 

r = ^< E ' I ^ II ^ I cos(0 4 -e j )>-k<[J2 1 1 h II to I MOi - o,)] 2 >}, (7) 
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where the primes denote that the sums are carried out over links along one of the 2 directions (x or y). 
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FIG. 1. (a): v vs. T for A = 0.01 (triangles up), A = 0.1 (triangles down), A = 1 (*), A = 2.5 (x), A = 10 
(+) and XY (circles), (b): zoom of l-(a) to compare the drop to of v for A = 0.01 with the less sharp drop for 
A = 0.1. 



Fig. T(a) shows a plot of v vs. T for different values of A and L = 40. For A = 0.01 (* = 0.94 we 
observe a sharp jump in the vortex density v (triangles up). As long as we increase A the jump becomes 
more smooth and moves to higher values of T c until for A = 10 (+ symbols) we get something very close 
to the A — > oo KT behavior (circles). The increase in the density of vortices when amplitude fluctuations 
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are large is in agreement with the analytical computations of ref. [|13|. What is new is the fact that, 
when £ ~ a (A = 0.01), the transition occurs almost directly from vortex to a plasma of vortices i.e. 
the bound pairs of vortex-antivortex seem to play no major role. Fig. l-(b) is a zoom of l-(a) showing 
the difference between A = 0.01 and A = 0.1 below T c : for A = 0.1 v drops to non-zero values due to 
the existence of vortex-antivortex pairs while for A = 0.01 v drops much more sharply to signaling the 
sudden extinction of vortices. 

The different behavior of v above T c between the large fluctuating amplitude LGW-regime - ~ 1 and 

the KT-regime * -C 1 is because amplitude fluctuations -governed by the A parameter- indeed decrease 
the energy of vortices enhancing vortex production. The same happens for the XY model with modified 
interaction jHj]; in fact, the shape modification can be straightforwardly connected to a core energy 
variation. 

The scarcity of bound pairs of vortex-antivortex below Tc for the extreme fluctuations regime (very 
small values of A or £ ~ a) can be explained in terms of the behavior of their free energy F pa i r = 
Epair — TSp a i r at T ~ T c from below. Roughly, E pa i r ~ 2E C , where E c is the vortex core energy, 
and Sp a i r ~ ln(-|j-). The three quantities S pa i r ,E c and T c all decrease as A decreases making difficult 
to disentangle the "energetic" contribution to v, proportional to exp[—E pa i r /T c ], from the "entropic" 
contribution ~ l/£ 2 . For the intermediate range of A (or £), although it is not easy to predict the 
behavior of the energetic factor, the entropy seems to be the main responsible for lowering the density 
of bound pairs |fl8|| . From a A sufficiently small, T c starts to decrease with A faster than E pa i r and thus 
exp[— Fp a i r /T] decreases more and more sharply making smaller and smaller the probability of bound 
pairs of vortex-antivortex. 

Figures 2-4 show plots of e, v and T for A = 0.01,0.1 and 10. For A = 0.01 the transition is clearly 
first-order: we observe latent heat and discontinuous changes in v and T at T — T c . On the other hand, 
for A = 10 (or | ~ ) the results are similar to those of the XY model. In particular, we get something 

close to the KT universal jump Ar = — . As long as A decreases we get a larger non universal jump in 
r. A = 0.1 corresponds to something in between first order and KT. 
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FIG. 2. e, v and r vs. T , A~T= 0.01 for sizes: L = 10 (circles), L = 20 (diamonds) and L = 40 (squares). Error 
bars for e and t> are smaller than the symbol sizes. 




FIG. 3. e, v and T vs. T, A = 0.1 for sizes: L = 10 (circles), L 
bars for e and w are smaller than the symbol sizes. 



20 (diamonds) and L = 40 (squares). Error 
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FIG. 4. e, w and Y vs.' T, A = 10 for sizes: L — 10 (circles) , L = 20 (diamonds) and L — 40 (squares) . Error 
bars for e and v are smaller than the symbol sizes. 



The double peak structure corresponding to the 2 coexisting phases, characteristic of a first-order tran- 
sition, is showed in Fig. 5(a) for A = 0.01 and sizes L — 10, 20 and 40. The width w of each of the peaks 
clearly scale asw~ \flj^n) = j; due to ordinary non-critical fluctuations. Fig 5(b) is a zoom of the right 
peak centered around the higher energy phase. It shows that w ~ 0.1(L = 40), w ~ 0.2(L = 20) and 
w ~ 0.4(L = 10). 

Therefore, one has a simple model in which the nature of the phase transition depends on the value 
of one parameter (A or £/a) which controls the thermal fluctuations of vortex cores and from which the 
following picture emerges: 

1) For A < 0.1 or | ~ 1 the density of vortices experiments an abrupt jump which coincides with a 
first order transition with large latent heat and a non universal jump in T all at a T c which grows with A. 
This is the LGW-regime as opposed to the much more smooth KT-regime. In particular, for sufficiently 
small values of A (for instance A = 0.01) below T c v drops to zero i.e. the number of measured bound 
pairs of vortex-antivortex is negligible compared with the number found in the KT transition. 

2) For A > 10 or (|) 2 < 1 we get basically the XY model (| ip | is fixed = 1 for all K) and the more 
subtle KT transition (with an unobservable essential singularity in the specific heat at T c and a much 
more small non-universal maximum above T c and a universal jump in V). The number of vortices and 
the energy evolve smoothly across the transition. 

3) For intermediate values of A (0.1 < A < 10) we have an interpolating regime between LGW and 

KT. 
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FIG. 5. (a) Histograms of e for A = 0.01 and L = 10, 20 and 40 at T = T c . (b) Zoom of the right peak. 



Whether or not a large enough increment of j or A to alter the nature of the phase transition driven 
by vortices can be accomplished by varying some thermodynamic parameter, for instance the pressure, 
is something which deserves investigation. 

This change in the nature of the phase transition when amplitude fluctuations of ip are not negligible 
(A < 0.1) is in agreement with very recent variational computations for the same model |l9| . Furthermore, 
we checked that all the couples of values (K c , A) for what we found a first order transition are such that 
equation (^) gives 0.7 < £/a < 1.1 in complete accordance with figure 6 of ref. (T^| from which we can 
see that for £/a bigger than 0.5 the RG trajectories cross the first order line of Minnhagen's |nj generic 
phase diagram for the two-dimensional Coulomb gas ]20|] . The jump in T larger than the universal value 
seems consistent with experiments on thin films HT C superconductors pl| . 

Finally, our analysis could shed light on the nature of the melting transition for the 2D FLL which 
remains controversial theoretically as well as experimentally. Experimental works and numerical simula- 
tions favor a KT-like transition in some cases and a discontinuous one in others p2"[ . A recent extensive 
Monte Carlo simulation found a first order transition at a temperature close to the estimated one assum- 
ing a KT melting transition ^3|. After all, it is possible that the nature of the melting transition in 2D 
strongly depends on the particular conditions and details of the studied specimen which in turn translate 
into different values of -. 

a 

Work supported in part by CSIC, Project No. 052 and PEDECIBA. We are indebted with D. Ariosa, 
H. Beck, P. Curty, E. Dagotto and G. Gonzalez-Sprinberg for valuable discussions. 



7 



[1] G. Blatter, M.V. Feigelman, V.B. Geshkenbein, A.I. Larkin and D. Vinokur, Rev. Mod. Phys. 66, 1125 (1994) 

and references therein. 
[2] V.L. Berezinskii, Sov. Phys. JETP 34, 610 (1972). 
[3] J.M. Kosterlitz and D.J. Thouless, J. Phys. C6, 1181 (1973) 
[4] M. R. Beasley, J.E. Mooij and T.P Orlando, Phys. Rev. Lett. 42, 1165 (1979). 
[5] S. Doniach and B.A. Huberman, Phys. Rev. Lett. 42 , 1169 (1979). 
[6] A.F. Hebard and A.T. Fiory, Phys. Rev. Lett. 50 , 1603 (1983). 

[7] J. Corson, R. Mallozzi, J. Orenstein, J.N. Eckstein and I. Bozovic, Nature 398, 221 (1999). 

[8] A. Jonsson, P. Minnhagen and M. Nylen, Phys. Rev. Lett. 70, 1327 (1993). 

[9] E. Domany, M. Schick and R. Swendsen, Phys. Rev. Lett. 52, 1535 (1984). 
[10] J. E. van Himbergen, Phys. Rev. Lett. 53, 5 (1984). 
[11] P. Minnagen and M. Wallin, Phys. Rev. B36, 5620 (1987). 
[12] S. E. Korshunov, Phys. Rev. B46, 6615 (1992). 
[13] D. Bormann and H. Beck, Jour. Stat. Phys. 76, 361 (1994). 
[14] M. N. Barber, Phys. Rep. 59, 375 (1980). 

[15] In fact this is the usual parameterization in lattice quantum field theory for the Euclidean action; see for 

instance K. Jansen, J. Jersak, C.B. Lang, T. Neuhaus and G. Vones, Nuc. Phys. B 265, 129 (1986). 
[16] M.E. Fisher, M.N. Barber and D. Jasnow, Phys. Rev. A8, 1111 (1973). 

[17] An analogous expression for the case of a granular superconductor appears in C. Ebner and D. Stroud, Phys. 

Rev. B28, 5053 (1983). In a similar way is calculated the orbital magnetic susceptibility, see for instance: A. 

Sewer, H. Beck, X. Zotos; Physica C 317-318 , 475 (1999). 
[18] We measured both E c and v for different values of £ at T ~ T c from below and found that the entropic factor 

fits in well with v in the intermediate range of £ while the energetic factor is the dominant one when £ ~ 1 

and v becomes neglegible. 
[19] H. Beck and P. Curty accepted for publication in Physical Review Letters. 

[20] Formulas (87) and (88) of |13| , connecting (£,T) with (Tcg,zcg), provide a link between our calculations and 
Minnhagen phase diagram. However, one should not compare directly values of (A, if) to (Tea, zcg), since 
the former are initial model parameters and the latter vary during renormalization according to the curves 
shown in B-B figure 6. 

[21] Ch. Leeman et al, Phys. Rev. Lett. 64, 3082 (1990). 

[22] K.J. Strandburg, Rev. Mod. Phys. 60, 161 (1988). 

[23] Y. Kato and N. Nagaosa, Phys. Rev. B47, 2932 (1993). 



8 



